Probability density of determinants of random matrices. 
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Abstract. In this brief note the probability density of a random real, com- 
plex and quaternion determinant is rederived using the singular values. The 
behaviour of suitably rescaled random determinants is studied in the limit of 
infinite order of the matrices. 



1 Introduction and results. 

We consider nxn matrices whose elements are either real, or complex or quater- 
nions (in what follows, the four components of the quaternions will always be 
real); the real parameters entering these elements are independent gaussian 
random variables with mean zero and the same variance. The number of real 
parameters needed to characterize an n x n matrix is thus f3n 2 , where (3 is 1, 

2 or 4 according as the matrix elements are real, complex or quaternions. We 
will derive the probability density of their determinant. 

The determinant of random real matrices is an old subject Q], that of ran- 
dom complex matrices and of random hermitian complex matrices was studied 
some time back ^ , that of random quaternion matrices presents some peculiar 
features due to the non-commutative multiplication, as we will sec below while 
the case of real symmetric matrices has been settled recently for odd n ||. 

The method we will use here is to start with the joint probability density of 
the singular values rather than that of the eigenvalues. As the absolute value 
of the determinant is the product of all the singular values, we can find its 
probability density by calculating its Mellin transform. This gives new proofs of 
the known results for random real and complex matrices and of a partial result 
for quaternion matrices. 
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For quaternions, multiplication being not commutative, it is not possible to 
define a determinant having the usual three properties Q ; namely, (i) det A = 
if and only if Ax — has a non-zero solution x ^ 0, (ii) det(AB) — det A. det B, 
(iii) det A is multi-linear in the rows of A. So the definition of a determinant 
varies according to which of the property or properties one wants to keep. We 
will adopt the following definition due to Dieudonne or Artin ||. 

Any matrix A is either singular (i.e. Ax = yA = have non-zero solutions) 
or has an inverse (i.e. AB — BA = I) 0. If A is singular, define det A = 0. If 
A has an inverse, define det A by recurrence on n as follows. If n = 1, define 
det A = | an |, where \x\ means the norm of (the quaternion) x. If n > 1, then 
let Aij be the (n — 1) X (n — 1) matrix obtained by removing the i-th row and 
the j-th column of A. The matrix elements of B, the inverse of A, are written 
as bij. Not all bij are zero. One shows || that whenever bij = 0, det A,-.; = 0, 
and whenever bij ^ 0, det Aji ^ and |&^~. det Aj%\ is independent of i or 
j. One then defines det A — \b^ det A,i|. Thus for a quaternion matrix A, 
det /I is a non-negative real number. (For real or complex A this definition 
also gives a non-negative real number, the absolute value of the usual ordinary 
determinant.) Note that this determinant is not linear in the rows of A, but has 
the other two properties |5|. Also that a quaternion matrix A may be singular 
while its transpose has an inverse [Q . 

The eigenvalues and eigenvectors of a matrix A are defined as the solutions 
of Aip — tpx, where ip is an n x 1 matrix and a; is a number. For a real or complex 
A one can eliminate ip to get det (A — xl) — 0, where I is the unit matrix. For 
quaternion A, if x is an eigenvalue with the eigenvector if and n any constant 
quaternion, then [i~ x x[i is an eigenvalue with the eigenvector ip/j,. Thus x and 
fi~ 1 xfi are not essentially distinct as eigenvalues. It is not evident that an n x n 
quaternion matrix should have n (quaternion) eigenvalues, but it has ^j. One 
can actually put them in correspondance with complex numbers 0. Here we 
will only note that the norm of the product of eigenvalues gives the determinant 
defined above. 

If all the eigenvalues of A are essentially distinct, then one can diagonalise 
A by a non-singular matrix. To make things clearer, we give an example: 
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showing that the eigenvalues of 



1 e 2 
ei e 3 



are and 1 — e%, while those of its 



transpose 



1 ex 

e2 e 3 



are x\ and X2- Their determinants are respectively and 2. 



If all the eigenvalues real and positive (respectively, real and non- 

negative), one says that A is positive definite (respectively, positive semi-definite). 
Denote by A* the transpose, hermitian conjugate or the dual of A according as 
A is real, complex or quaternion. For any matrix A, the product AA* (or A^A) 
is positive semi-definite, its eigenvalues are real and non-negative. The positive 
square roots of the eigenvalues of AA* (or of A^A, they are the same) are known 
as the singular values |j| of A. The eigenvalues and singular values of A have, 
in general, nothing in common, except that 



Jl = det(AA t ) = det A. det A^ = | det A\ 2 = f[ 



(1.5) 



where Xi are the singular values and the eigenvalues of A. 

In section 2 we start with the joint probability density of the singular val- 
ues and calculate the Mellin transform of the probability density p(|y|) of the 
(absolute value) of the determinant y of a random matrix A. From symmetry, 
when A is real, y is real and p(y) is even in y\ when A is complex, y is complex 
and p(y) depends only on \y\. When A is quaternion, y is by definition real 
and non-negative. One can therefore recover p{y) from p(|y|) when A is real or 
complex. Our results, confirming those in the known cases, are as follows. 



pi(y) = IlFC?/ 2 )]" 1 ^ 
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(1.6) 



(1.7) 
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,2n 

2 



y real non negative. 

(1.8) 



Here Gjjf is a Meijer G-function. In the above results the gaussian proba- 
bility distribution P(A) for the matrix A was taken P(A) cx e" atrA A , with 
a = 1. Next we show that the probability density of the random variable y = 
[det A'l'A] 1 /" converges to 5{y— 1/e) in the large n limit , with P(A) cx e ~ atrA ^ A 
and a = (3 n/2. In section 3 we study the large n limit for a non-gaussian random 
complex matrix and show that the random variable y — [det A^ A] 1 / 71 converges 
in the large n limit to a constant, whose value, depending on the parameters 
in the non-gaussian probability distribution, is different in the two phases of 



3 



the model. Finally for a hermitian complex random matrix H with probability 
density P{H) cx e ~ ntrH , we show that in the large n limit | det 7J| 1//n tends to 
the constant l/v2e- 

Some of these results are probably known to some experts, since analogous 
results appear in the literature ||. 



2 Gaussian matrices. 

The joint probability density of the singular values can conveniently be derived 
in two steps from the two observations (l(J , 

(i) Any matrix A can almost uniquely be written as UAV where A is a 
diagonal matrix with real non-negative diagonal elements, while U and V are 
real orthogonal, complex unitary or quaternion symplectic matrices according 
as A is real, complex or quaternion; "almost uniquely" refcring to the fact that 
either U or V is undetermined up to multiplication by a diagonal matrix. 

(ii) Any positive semi-definite matrix H = AA^ can be written uniquely 
as H = TT\ where T is a triangular matrix with real non-negative diagonal 
elements. 

As a result the gaussian joint probability density exp(— atiAA^) for the 
matrix elements of A gets transformed to 

F(A) ee F(\i, A„) = const, exp | -a^Aj J | A(A 2 )|^ f[ Aj ^ (2.1) 

where Ai, .., A„ are the singular values of A, A is the product of differences 

A(A 2 )= J] (\l-X]) (2.2) 

l<j<k<n 

and /3 = 1, 2 or 4 according as A is real, complex or quaternion. 
The Mellin transform of the product of the A's is 



M n (s) = const. / r] s S(rj - Ai...A„) F(A) dXi...d\ n dr\ 
Jo 

r oc I n \ ™ 

= const. J exp -aj^ A 2 |A(A 2 )|^ J] X^+^dX, 



3=1 / 3=1 



r oo ( n \ n 

const. J exp | -aj^ *i j ^(i)!' 3 J] tf + "~ 3)/2 #; 
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3=1 

In the last line we have used a result derived from Selberg's integral |t^]. The 
constant, independent of s, is fixed from the requirement that M n (l) = 1. 

The inverse Mellin transform of the expression (^^) is a Meijer G-function 



I s-l i j/3 
\ 2 ' 2 



r(f 



(2.3) 



P p(\y\) = 2a"/ 2 J] [r( J /3/2)]- 1 G$ ( a»|y| 
i=i 



2/3-1 n/3-T 



(2.4) 



When /3 = 1, the matrix A is real, its determinant y is real, from symmetry the 
probability density p\ (y) is an even function of y and we have 



pi(y) = 2^i(bl) 



(2.5) 



giving equation (1.6) with a — 1 . When /? = 2, A is complex, ?/ is complex, 
from symmetry P2(y) depends only on the absolute value \y\ of y, and one has 
for a = 1 , 



P2(y) = 



l 



-M\v\) 



3 = 1 
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,n - 1) 



(2.6) 



3=1 



which is equation (1.7). When /3 = 4, v4 is quaternion, y is, by definition, real 
positive, and Pi{y) — Pi{\y\), giving equation (1.8). 

In reference 2, appendix A. 5, we somewhat conventionally mapped the (quater- 
nion) eigenvalues on to the essentially equal eigenvalues having the scalar part 
and only one other component at most, accounting for a factor y 2 in the proba- 
bility density. Moreover, equation (A. 43) there has a misprint, T((s + 2j + l)/2) 
there should read T((s/2) + 2j + 1). Thus equation (1.8) tallies with equation 
(A. 43) of reference 2. 

We now evaluate the large n behaviour of the moments < y k > of the 
random variable y — [dei A 1 A] 1 /™ which show that it converges to a constant in 
the n — > oo limit. Next the same result is obtained by the saddle point method. 
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In the study of large n, the proper choice of the parameter a for the gaussian 
ensembles is a — (3n/2. Then eq.(2.3) implies 



< detAU >= 



(3n 



zk n 



n 

i=i 



r(fc + f ) 
r(f) 



(2.7) 



that is 



( \ 3n T(k + i£) 

log < det A^A >= -nk log -— + log - 2 



rtfc log — — h « / ax log 



r(f) 

h 

2 ; + 0(log n) 



/3 _l_ 7ix/3 ^ 



> 2 ' 2~ 



—nk + 0(log n) 



(2.8) 



where the Euler-Maclaurin formula has been used to estimate the large n asymp- 
totics. When x is near 0, the integrand is a constant and its contribution is neg- 
ligible. When x is not small, o ne c a n ig nore other terms compared to nx(3/2. 
Replacing k — > k/n in eqs.(2/7), (|2.8|), one gets 



k/n 



lim log < det A*A >= -k 



(2.9) 



From the knowledge of all the moments (2.9), we conclude that the random 
variable y — (det A f A) 1 ' n converges in the large n limit to the constant 1/e. 



It is convenient to evaluate the above large n limit also by a saddle point 
approximation because this is easy to generalize to different probability distri- 
butions. Let us recall the asymptotic density of squared singular values (see |Q 
eq.(10) after setting L = 1, m 2 = 1, g = 0, hence A = 0, B = 4 ) 



pit) = lim - V < 5(t - tj) >- 



1 A-t 
2^V t 



< t < 4 



(2.10) 



Then it is easy to evaluate 



lim 

n — >oo 



1 ( ^ 1 

- < log detA^A >= lim - < /c^log t, > 
n \ / n— >oo n z — ' 

\ / j=i 



\ogtp(t)dt 



(2.11) 



confirming eq.rt2 



G 



3 Large n for non-gaussian complex matrices. 

Let us now consider an example of non-gaussian probability distribution such 
that in the large n limit two different spectral densities for the singular values 
exist. For simplicity, we consider the ensemble ofnxu complex matrices A 
with the non-gaussian probability distribution 

P(A) cx e -n(otrAtA+26trAtAA+A) f b>0 , a real (3.1) 

The analogous ensemble of real matrices would need only trivial changes. 

Again the evaluation of all the moments of y = [det A^ A] 1 ^ 71 may be per- 
formed in terms of tj , the squared singular values of the matrix A 

Since Y[j=A^j) k ^ n — ex P 1°£^) > ^he large n limit for all the moments 

< y k > are easily evaluated by the saddle point approximation 




lim < y k > = 



lim < exp ^— ^logt,^ >= exp J dt p(t) log tj (3.3) 



where p(t) is the solution of the saddle point equation 

a 



Z +2 bt= i^-dy (3.4) 
2 J t-y 



The solution of eq.(3.4) has two different forms, pi(t) , p2(t) (see eq.(10) 
after setting L = 1, m 2 — a, g — b ), depending on the values of the real number 
a/Vb being larger or smaller than the critical value a/Vb = —4. At the critical 
value, px(t) = p 2 {t). 



M 1 C — t 

7T V t 
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< t < C 



Va 2 + 486 - a r- 
C = — , a/Vb>-4 ; (3.5) 

The definite integrals related to eq. (|]^) for the spectral function pi (t) are known 
in closed form and 



c 



C 1 ciC 

pi (t)logtdt = log----— (3.6) 



Therefore in the "perturbative phase", that is for a/Vb > —4 , the random 
variable y converges in the large n limit to a constant 

lim [det AU] 1 /" = - e-i-^r (3.7) 

n — >oo 4 
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For a/yb < —4 we have 



2 b 



P2 (t) = -y/(B-t)(t-A) 

7T 



A < t < B 



A + B = - 



2b 



B-A=— , a/Vb<-4 
Vb 



(3.8) 



The definite integral related to eq^^) for the spectral function p 2 {t) may still 
be evaluated 



I(a,b) 



p 2 (t) logtdt 



2b 



I" J A 



yJ{B -i)(t- A) log t dt 



= b -(VB-VAr+ b -(B-Aflog 



A 



(3.9) 



One may still conclude that also in the "non-perturbative phase" , that is for 
a/Vb < —4 , the random variable y converges in the large n limit to a constant 



lim [det^ f ^l] 



l/n _ J{a,b) 



i/Vb < -4 



(3.10) 



We remark that other functions of the determinant, such as w — [det A^yl]", 
may not have a finite limiting probability distribution, in the large n limit, 
though their large n behaviour may still be evaluated. 



Hermitian matrices. The large n limit of the absolute value of the 
determinant of gaussian hermitian matrices may be evaluated from the exact 
finite n moments given in [Q. Let us consider an ensemble of n x n hermitian 
matrices H with probability distribution P(H) oc e ~ ntrH and the random 
variable y 



V 



I det H\ 1/n 



Let Xj , j = 1, .., n be the eigenvalues of the matrix H, then 

k /..r 00 ^(,)n; =1 Kf/ n e-«4d, j 

<y >= 



These moments are known from Si. 



n p/l 



2n 



(3.11) 



(3.12) 



(3.13) 



log < y k >r 



-|log( 2e ) + (^ 
2 \ n 



(3.14) 
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Thus in the large n limit | det H\ x / n tends to the constant l/y/2e. 
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